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Nonequilibrium electron distribution in a superlattice subjected to a homogeneous electric field 
^\ . (biased superlattice with equipopulated levels) is studied within the tight-binding approximation, 

taking into account the scattering by optical and acoustic phonons and by lateral disorder. It 
is found that the distribution versus the in-plane kinetic energy depends essentially on the ratio 
between the Bloch energy, eb, and the optical phonon energy, hu>o- The in-plane conductivity is 
calculated for low-doped structures at temperatures 4.2 K and 20 K. The negative conductivity is 
found for bias voltages corresponding to es/fiojo —1/2, 1/3, 2/3... (the Bloch-phonon resonance 
condition) . 
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I. INTRODUCTION 



Vertical charge transfer in a superlattice subjected to a homogeneous electric field (biased superlattice, BSL, with 
equipopulated levels) has been under investigation starting 70th (see Ref. [l[ and references in the reviews of Ref. [H). 
The stimulated emission in the mid-infrared (IR) and terahertz (THz) spectral regions, caused by the intersubband 
transitions of electrons under vertical transport through tunnel-coupled cascade structures (monopolar laser effect), 
has also been investigated. Using this scheme, both mid-IR and THz lasers viability has been demonstrated during 
the previous decade (see Refs.fUfJ] and references therein ). Recently, nonequilibrium electron distribution has been 
observed experimentally [f| and described theoretically @ for heavily-doped cascade structures, when the effective 
temperature is determined from the balance equation. To the best of our knowledge, there is no consideration of 
1 ^ 1 ' nonequilibrium carriers for low-doped structures which were performed beyond the balance approach. 

In this paper we study the nonequilibrium electron distribution in a biased superlattice (BSL) under vertical current 
through the Wannier-Stark ladder, which takes place under the condition 2T <C eb, where Eb is the Bloch energy 
and T stands for the tunneling matrix element between adjacent quantum wells (QWs) Q. Since the parameters of 
each QW and the conditions for interwell transitions are identical (see Fig. la) the level populations over QWs are 

^ I the same. But the distribution of electrons over the in-plane energy should change essentially due to the interplay 

. between elastic and non-elastic processes (see Figs. lb-d). In the low-temperature case we only consider the passive 
region, with energy less than the optical phonon energy, Kloq. For the low-concentration limit, we consider the kinetic 
equation which takes into account the quasi-elastic scattering caused by acoustic phonons as well as the interwell 
tunneling due to elastic scattering by disorder and due to optical phonon emission. As a result, we obtain the electron 
distribution versus the in-plane kinetic energy which strongly depends on the ratio £s/?ia;o- In the case of the Bloch- 
phonon resonance, when Meb = NHluq with N and M integers, a partially-inverted distribution, with maxima at 
. £h ! energies (N/M)Hujo, can be realized. 

The phenomenon of absolute negative conductivity (ANC) of electrons excited near the energy hivo was discussed 
four decades ago [|[ and different regimes of ANC (including magnetotransportQ and transient regimes of the 
response (To|) were considered. Recently, ANC reg ime was observed when microwave radiation acts on two-dimensional 
(2D) electrons in a quantizing magnetic field [11]. As shown below, a resonant ANC regime of the in-plane response 
can be obtained in BSL under the Bloch-phonon resonance conditions. Such a peculiarity appears due to the contri- 
bution of the energy interval near huo where an inverted distribution takes place. As a result, BSL becomes instable 
with respect to in-plane fluctuations if ANC conditions are satisfied. 

The paper is organized as follows. The basic equations describing distribution of hot electrons in BSL and in-plane 
conductivity are presented in Sec. II. Analytical consideration for the case Eb/^o = 1/2 (the second-order Bloch- 
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FIG. 1: (a) Band diagram for BSL of period Z with the Bloch energy eb which is comparable to the optical phonon energy, 
Hljo. (b-d) Schemes of tunneling transitions due to elastic scattering (solid arrows), spontaneous optical phonon emission 
(dotted arrow), and phonon emission from the active region (vertical dashed arrows) for the cases: (b) eb < fkJo/2 , (c) 
< eb < hivo , and (d) Hloq < eb- 

phonon resonance) is performed in Sec. III. Results of numerical calculations are discussed in Sec. IV. Concluding 
remarks and the list of assumptions made are given in the last section. In Appendix, the kinetic equations for different 
es/frujo are presented. 

II. BASIC EQUATIONS 

Within the tight-binding approximation, the electrons in BSL are characterized by the 2D momentum, p, and the 
quantum well number, r = 0, ±1, . . .. Under the in-plane electric field E, the distribution function, f rp , is governed 
by the system of kinetic equations: 

k 

where the collision integrals J/-(/|rp) describe the scattering processes caused by the longitudinal optical phonon 
emission (k — LO), the acoustic phonons (k = ac), or the static disorder (k = d). Below we present these collision 
integrals and consider the kinetic equation for the distribution functions f rp , which are normalized by the condition 
n = (2/V) J2rp frp: where n is the 3D concentration, V is the normalization volume, and the factor 2 is due to spin. 
We also evaluate the lateral current density I = (2e/Vm) J2 rp P/rp i0r electrons, with the effective mass m, under a 
weak probe field E. 

A. Collision integrals 

Here we evaluate the collision integrals in Eq. (1) by modifying the general expressions [l2| for electrons in BSL, 
described in the tight-binding approximation by the states \rp), with the energies s r p — v*&b £p: where Sp — f^vn 
is the in-plane kinetic energy (see Sec. 3 in Ref. 7). For the low-temperature case, when temperature of phonons 
T p h <C hu>o, the spontaneous emission of dispersionless optical phonons is described by 

JLo(f\rp) = y \C { Q LO) \ 2 \(r'p'\e lQr \rp)\ 2 [5{e r , p , - e rp - huj )f r>pr - S(e rp - £ rV - Tiuj )f rp ] , (2) 

r'p'Q 
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where ICq^ ^ 2 is the bulk matrix element for the Frolich interaction, with the vibration mode characterized by the 
3D wave vector Q, and |(r'p'|e lQr |rp)| 2 is the overlap factor. Taking into account the quasielastic energy relaxation 
caused by the equipopulated acoustic phonons, one obtains the collision integral 



Jac(f\rp) = Wrpr'p'(/r'p' - f rp ) - ]- ^ A W rprV (f r , p . + f rp ). 



(3) 



r'p' 



r'p' 



The transition probabilities W rpr ' p > and AW rpr / p / are written here within the second order accuracy with respect to 
the acoustic phonon energy, Hcoq, as follows 



K rr ,{p - p')5(e r > p i - e rp ) + -^-AK rr >(p - p')5"(e 



r'p' ^rpj-, 



AW rpr > p > = AK rr ,(p - p')5(e r > p > - s rp ). 



(4) 



Here S'(E) and 5"{E) are the first and second derivatives of the 8- function and the kernels K rr > and AK rr ' are given 

by 



K r AP-p') = ^£l4 ac) |VpV Q Hn>)l 2 |f 

Q 



Q 



AK rr ,(p-p') = ^E|cJ c) | 2 |(r'p'| e ^|^P)| 2 ^ Q 



(5) 



Q 

where Cq°^ is the bulk matrix element for the deformation interaction. 

We restrict ourselves to the sequential tunneling processes under the condition T <C £b- Considering only the 
proportional to (T/eb) 2 corrections to the overlap factors, we use 



T 

EB 



(Srr' + l + ^rr'-l) 



(6) 



where 1 & q± d = \(0\e tq±z \0)\ 2 describes the transverse overlap between the ground states of the QWs, |0). Since all 
QWs are identical and s r i p i — £ rp = (r — t')eb + E p > — E p , the distribution functions are the same in any QW, i.e. 
frp — * f P - Thus, the collision integrals in Eqs. (2) and (3) are independent on r because the summation over r' is 
replaced by X^Ar=±i ■ ■ ■• The collision integral in Eq. (2) is transformed into 



r'p'q± 



X <^ 5{£p> - £ p - ?lU>o)fp' - S(e p - E p ' - hu> )f p 



+ 



Ar=±l 

-S(Ar£ B + 



5{Ar£ B + Ep' - £ p - hio )f p ' 

}■ 



(7) 



where X^Ar=±i ■ ■ ■ describes the interwell tunneling with LO-phonon emission and Q 2 = \p — p'| 2 /?i 2 + q\. 

Below we restrict ourselves to the thin QW case, when |Cq C ^| 2 can be replaced by |Cg"°^| 2 . Similar transformations 
for the acoustic phonon contribution of Eq. (3) give us 



Jac{f\p) = K ac^2 H £ p'- e p) + E [j^l S ( Ar£ B + Ep' - S p ) 



Ar=±l 



lfp> fp) 



-AKJ2^S"(e p , E p ){f p , / p ) - 1 -5'(E P ' E p )(f p , + f p ). 



(8) 



p' 
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Here we have neglected weak (oc AK) contributions to the tunneling transitions. The kernels in Eq. (5) appear to be 
momentum independent 



4tt 



T, 



ph 



Two 



q± 



AK « ^El^ c) l%x^ 



(9) 



q± 



due to the narrow QW approximation. The intra- and inter-well scattering caused by the static disorder can be 
described in a similar way to the elastic (oc K ac ) contributions in Eq. (8): 



Jd{f\p) = K d^2 5 ( £ p'~ £ p) + ^2 ( — ) 5{Are B + e p , - e p ) 

p' <- p'Ar=±l \ £B 



Up' - / P ) 



(10) 



Factors Kd and K ac determine the departure relaxation rates caused by the elastic scattering mechanisms as Vd,, 
Kd,acJ2 P ' fi( £ p' ~ %>) P2D1 where p 2D is the 2D density of states. 



B. Nonequilibrium distribution 

We search for the solution of Eq. (1) in the form / p ~ f £ + Af p , where f e describes the lateral heating due to 
tunneling current and A/ p is the in-plane anisotropic addendum due to the weak field E. We consider the symmetric 
part of the distribution which is governed by the kinetic equation J2k Jk{f\s) — and satisfies the normalization 
condition 

/>OC 

nZ = p 2D / def e (11) 
Jo 

with the layer concentration nZ. Averaging Eq. (7) over p-plane and taking into account the energy conservation 
condition, one obtains the LO-contribution as a finite-difference form: 

JLo(f\z) — V £ +hujo,efe+hu>o ~ v e-huj ,efe + ^ l^e+Swo-ArEB ^fe+huj a -Are B ~ -hu> - Are B ,efe] ■ (12) 

Ar=±l 

Here the tunneling contributions v E = {T/£b) 2 ve,c are reduced by the factor (T/eb) 2 - We introduce the relaxation 
rate describing the spontaneous emission of LO-phonons as follows 

ve, e = y E \C { Q LO) \ 2 * q±d 6(e p ,~E). (13) 
p'ix 

Performing the integration over p'-plane one obtains 



" e(E)< ™° / **■ V (^TTTW^e - (14) 

— OG 

where a is the polaron coupling constant and Ed — (h/d) 2 /2m. This rate is of the order of auj . Fig. 2 shows the 
dimensionless relaxation rate VE,e/ a ^o versus e/huj , plotted in the passive region for different E/TujJ values and for 
a 60 A wide QW, when e d /huj ~0.44. Note, that ve^ appears to be logarithmically divergent if E — ► e. 

The intrawell process of quasi-elastic energy relaxation is described by the Fokker-Planck collision integral [l2| given 

by 

= 2 d ( d fe , fe 



Mf\s)^^-^ + ^J, (15) 

where v ac is the above-introduced departure relaxation rate and e 2 ~ (AK / K ac )T p h/2. The elastic tunneling relax- 
ation caused by disorder and acoustic phonons [see Eqs. (8) and (10)] is governed by the finite-difference contribution 

Jt(f\e)~y t £ 0(e - AresXiUArea - /*) (16) 

Ar=±l 
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FIG. 2: (Color online) Dimensionless relaxation rate i/E,e/otu)o versus e/htoo plotted in the passive region for E/huo = 0, 0.2, 
0.4, 0.6, 0.8 and 1. 



with the tunneling rate v t = (T/eB) 2 (vd + v ac ). Thus, the distribution f £ is governed by the equation 

Jac(f\e) + Mf\s) + J LO (f\e) = 0. (17) 

Moreover, in the active region, e > Tvuj , the main contribution is due to the spontaneous emission of LO-phonons 
[first and second terms of Eq. (12)]. 

In the active region, e > hivo, the distribution decreases fast. Thus, in the kinetic equation, one has to take into 
account the second derivative from Eq. (15) and the spontaneous emission contribution form Eq. (12): 

V 2 0-^o/, = O, (18) 

where v^o — ^o,7iw - Using the boundary condition f £ ^oo — one obtains the solution f e « f% Uo exp[— Ao(e — Thoq)] 
with A = v ' v 'lo i 'v ac^ 1 ■ 

Next, eliminating the fast spontaneous emission of LO-phonons one obtains the Eq. (17) in the passive region with 
the boundary condition 

(JL + A °) f ^ R ^° = °- (i9) 

Thus, the problem formulated in the passive region takes into account the quasi-elastic energy relaxation described 
by Eq. (15) and the interwell tunneling transitions shown in Figs. lb-d. The normalization condition should be 

restricted over the passive region and Eq. (11) is transformed into nZ = p 2D J Q dsf e . 



C. Lateral conductivity 

Further, we turn to the description of the linear response given by A/ p oc E and consider the current density 

2e f dp 



ZmJ (2^ Af v { ^ 
The nonsymmetric part of the distribution function A/ p is determined from the linearized kinetic equation 
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where the elastic collision integrals due to etc- and (^-contributions can be replaced by —is m Af p . Here v m = v& + v ac 
is the momentum relaxation rate due to the elastic scattering [see Eqs. (8) and (10)]. The non-elastic momentum 
relaxation due to the optical-phonon-induced interwell transitions, Jlo(A/|p), is given by the oc (T/eb) 2 contribution 
of Eq. (7). Introducing the energy-dependent function x e according to A/ p = (e/m)(E • p)x eI we transform Eq. (21) 
into the finite-difference equation 

= —VmXe + Wt+hwo-AreB^Xe+fruo-AreB ~ u t-hu -Are B ,eXs\- (22) 

Ar=±l 

Here v E e — (T/eB) 2 VE,e is the tunneling-induced relaxation rate where 

?E > £ = Y J2 \ C Q G) I 2 *^ d cos(p75')^(e - e P > - E) (23) 
p'q± 

which uses a similar notation to Eq. (13). 

Introducing the in-plane conductivity, er, according to I = crE we obtain 



hu>o 

e2 P2D 
mZ 



dsEX £ - (24) 



Here we neglect the contribution from e > Tiluq because of the smallness of x e — v lo(~ df £ /de) in the active region. 
Under the condition v m 3> v f LO Eq. (22) gives x £ — v Vn{^df E /de) and the conductivity takes the form 

/ p 2D hu \ e 2 n 2D . . 

a = a 1 /n Wo , cr = -. (25) 

\ n 2 D J mv m Z 

As a result, er/ero i s expressed through the distribution function at e = ?io;o and a negative lateral conductivity takes 
place at f huo > n 2 D/p 2 D huj a- 



III. SECOND ORDER BLOCH-PHONON RESONANCE 



Before analyzing the general problem, we consider the simple resonant case 2eb = hwo, when the distribution can 
be considered over two intervals: < e < Tiujq/2 and Tiljq/2 < e < Hluq. Introducing the functions over the interval 
< e < hu /2 according to f\ e = f e and f 2e = f e +hu , we transform Eq. (17) [or Eqs. (Al) and (A3) in Appendix] 
into the system 

Jac{h\e) + Vt (2/3e - he) + ^ +huo /2,ef2e = (26) 
Jac{h\e) + Vt (fie ~ 2/ 2e ) ~ vl £+huJa/2 f 2£ = 

Here v t — (T '/ 'eB) 2 v m is the elastic tunneling rate and v E was introduced in Eq. (12). The two second-order 
differential equations (26) should be solved with the boundary condition of Eq. (19), the normalization condition 

nZ = p 2D $^ o/2 de(fu + f 2e ), as well as the inhomogeneous conditions f\nu /2 = f2e=o and {df lE /de) hu]o/2 = 
(df 2e /de) . 

Fig. 3 shows the distribution function f £ versus in-plane kinetic energy obtained from the system (26) for different 
T p h, T, and v m . As can be seen, peaks of f s widen as temperature increases, reducing their maxima. The effect of 
interwell coupling oc T and elastic scattering oc v m is also evident, founding that peaks increase as these parameters 
do. Calculations have been made for a concentration n 2 D = 10 9 cm" 2 (or n = 10 15 cm~ 3 if Z =100 A), so that 
fie < 0.1 and electrons are non-degenerate. [13| For this value, n 2 B / p 2D Tiluq ~ 10~ 3 . Peaks located at Tiujo/2 and 
hujQ are of the order of 10~ 2 and 10 -3 , respectively. Therefore, the last value is big enough to obtain ANC because of 
fhui > n 2 o/ p 2D hujQ, according to Eq. (25). In order to magnify peaks at e /hajQ=0.5 and 1, we have limited vertical 
axis size. 

Fig. 4 describes normalized lateral conductivity <t/<7o given by Eq. (25) vs momentum relaxation rate v m , for 
temperatures T p h =4.2 K and 20 K, and tunneling coupling T — 5 and 3.5 meV. Because of the distribution function 
peaks at Tuvq increase as T and v m do (see Fig. 3), lateral conductivity decreases correspondingly leading to negative 
values for a wide region of parameters. As we saw in Fig. 3, the effect of the temperature is opposed to the previous 
ones. For high temperatures a increases to reach ctq and the possibility of having ANC disappears. 
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FIG. 3: (Color online) Distribution function f s vs in-plane kinetic energy obtained from the system (26) at T p h =4.2 K and 
20 K (solid and dashed curves) for (a) T=5 meV and v m — l ps _1 , (b) T=5 meV and ^ m =0.5 ps _1 , (c) T=3.5 meV and v m =l 
ps _1 , and (d) T=3.5 meV and ^=0.5 ps _1 . 



IV. NUMERICAL RESULTS 



We turn now to the numerical calculations of the nonequilibrium distribution f e governed by the Eq. (17), with 
the boundary condition defined by Eq. (19), and the normalization requirement [the explicit form of Eq. (17) for the 
cases Sb < hu>o/2, Hujq/2 < eb < hioo, and huo < Eb are given in Appendix]. We also analyze the lateral conductivity 
solving the finite-difference Eq. (22) [see explicit expressions (A2), (A4), and (A6)] and performing the integration 
in Eq. (24). Calculations below are performed for the GaAs/Alo.3Gao.7As-based SL, formed by 60 A wide QWs 
separated by barriers of 32 A (or 37 A) wide, which correspond to the tunneling matrix element T— 5 meV (or 
3.5 meV). We consider temperatures of T p h =4.2 K and 20 K as well as the effect of the elastic scattering variation 
through the momentum relaxation rate v m =l ps _1 and 0.5 ps _1 . It is convenient to use n^D = 10 9 cm~ 2 in spite of 
ct/cto does not depend on concentration. 

Fig. 5 displays the distribution function f £ vs the in-plane kinetic energy for different Eb/^o values (1/3, 1/2, 
2/3, and 4/3) obtained from the general Eq. (17) (see also Appendix). When comparing the case es/?ia; =1/2 of 
Fig. 5(a,b) with the panel (a) in Fig. 3, calculated in the former section for the same parameters, a good agreement 
is found. As mentioned before, temperature effect is reflected as a widening and decreasing of the peaks. For other 
N/M values f s shows lower relative maxima. 

Next we calculate the normalized lateral conductivity solving Eq. (22) and using f E obtained before. Fig. 6 
represents a /oq as function of the Bloch energy sb/Ju^o and for different temperature, elastic scattering and tunneling 
coupling values. General behavior shows a pronounced relative minimum located at Eb/Jt-^q =1, followed by other 
relative minima at 1/2, 1/3, 2/3, 1/4,... (in decreasing order). In the active region, when N > M, these peaks 
are practically negligible. Depending on parameters some of the peaks reach negative values. To clarify the effect of 
these parameters we can compare in pairs the panels in Fig. 6. Thus, comparing panels (a) and (b) we can see the 
temperature effect, which is similar to the one found before for the distribution function: peaks are wider and less 
pronounced, leading to less negative values. Comparing panels (a) and (c) the effect of the tunneling coupling can be 
visible: if T decreases (increasing barriers width) size of peak minima are reduced. Finally, an analogous behavior is 
observed comparing panels (a) and (d) to see the elastic scattering effect. Reducing v m , we obtain a decreasing of the 
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FIG. 4: (Color online) Normalized conductivity <j/<j , given by Eq. (25), vs momentum relaxation rate v m , for different 
temperatures and tunneling coupling values (T v h and T). Solid line: (4.2 K and 5 meV). Dashed line: (20 K and 5 meV). 
Dotted line: (4.2 K and 3.5 meV). Dash-dotted line: (20 K and 3.5 meV). 

peaks in a similar way. One can conclude that the most favorable conditions to get negative conductivity correspond 
to low temperature, and big tunneling coupling and elastic scattering values. 

In order to detail the shape of the normalized conductivity peaks we present some of them in Fig. 7. We have 
chosen the more noticeable ones, corresponding to Bloch energy sb/Ju^o around 1/3 [panel (a)], 1/2 (b), 2/3 (c), and 
1 (d), for T =5 meV, with temperature and elastic scattering values corresponding to panels (a), (b) and (d) included 
in Fig. 6. A breakdown of the peak symmetry is observable when eg/fiuo increases going from a quite symmetric 
peak for eb/^ujq close to 1/3, to a clearly asymmetric peak around 1. 

V. CONCLUDING REMARKS 

In summary, we have demonstrated that the negative lateral conductivity regime is possible in low-doped biased 
superlattices under the Bloch-phonon resonance conditions. When analyzing the dependence of a vs bias voltage, 
narrow negative peaks take place if Eb/^oJo is close to the ratio 1/2, 1/3, 2/3. . .. ANC regime appears to be most 
pronounced for low-temperature region in BSL with an effective interwell coupling. 

Next, we list the assumptions used. The main restriction of the ANC regime consists in neglecting the electron- 
electron scattering, imposed by the Maxwell distribution, with an effective temperature suppressing the high-energy 
part of the distribution. The condition a < does not depend on concentration for non-degenerate electrons. Thus, we 
have used in the calculations n 2 D < 10 10 cm -2 , where the electron-electron scattering is not the main scattering process 
(see Ref. 14, where different systems were considered). The evaluation of the limiting concentration, which requires to 
involve the electron-electron collision integral in Eq. (17), lies beyond the scope of this paper. Another approximation 
we have made is rather standard. In order to estimate the coefficients of the kinetic equation, we have used a tight- 
binding approach for the description of the electronic states 0- The use of the bulk model for phonon dispersion and 
electron-phonon interaction is a reasonable approximation for the GaAs/AlGaAs-based structures [15]. We consider 
the momentum relaxation due to short-range scattering neglecting a large-scale potential; last contributions require 
a special attention in analogy with the case of a single low-doped well. We restrict ourselves to the case of uniform 
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FIG. 5: (Color online) Distribution function f e vs in-plane kinetic energy for different eb/Suo values (1/3, 1/2, 2/3, and 4/3) 
in GaAs/Alo.3Gao.7As-based BSL at different temperatures (a) T p h=4.2 K, and (b) T p h=20 K. 
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FIG. 6: Normalized lateral conductivity vs Bloch energy eb/Tiuio and for different temperature, elastic scattering and tunneling 
coupling (T ph , v m , and T). (a): (4.2 K, 1 ps" 1 , and 5 meV). (b): (20 K, 1 ps" 1 , and 5 meV). (c): (4.2 K, 1 ps~\ and 3.5 
meV). (d): (4.2 K, 0.5 ps" 1 , and 5 meV). 



bias fields and equipopulatcd wells neglecting the possible domain formation (one can avoid instabilities of vertical 
current in a short enough BSLflfjj]). 

We should also mention that an experimental task for measure the lateral conductivity of BSL is not simple because 
a complicate contribution of the corner regions is possible. But, instead of the dc current measurements, one can 
use a high-frequency contactless study of the response using a transverse capacitor. In addition, under the instability 
conditions (if a < 0) a direct measurement of lateral conductivity is not necessary because the vertical current appears 
to be unstable. A detailed description of this unstable response requires a special consideration. 

To conclude, a low-doped BSL at low temperature is a suitable structure in order to realize the absolute negative 
conductivity regime. In addition, a similar behavior is possible not only for the BSL under consideration but also for 
the more complicated tunnel-coupled structure used in quantum cascade lasers. An instability of such a structure for 
the case of low doping and temperature is possible and should be checked additionally. 



APPENDIX A: KINETIC EQUATIONS 



Below we rewrite Eqs. (17) and (22) for the cases Eb < hui /2 (A), hu n /2 < eb < (B), and Tiojq < Eb (C). 
For the case A, when eb < frw — £b < (see Fig. lb), the distribution in the passive region f e is described by 
the system 

Jac(f\s) + V t {f e +e B + fe+Hu -e B ~ fe) + ^l +hLJ0 _ £B £ fe+hij n -eB = °> < £ < £ B 

Jac(f\£) + Mfe+es + fe+-e B ~ 2 fs) = 0, £ B < £ < TiUJq - E B (Al) 

Jac(f\£) + Mfe-es ~ 2 fe) ~ ^t-Hu +e B ,efe = °) kuJ - £ B < £ < TlL) 

which is written for the three intervals < s < £b, £b < £ < ThjJq ~ e B, and Hluq — Eb < £ < hajQ. Eq.(22) in this 
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FIG. 7: (Color online) Normalized conductivity vs Bloch energy eb/Tt<uo around 1/3 (a), 1/2 (b), 2/3 (c), and 1 (d), for T =5 
meV. Solid lines: u m =l ps and T p h =4.2 K. Dashed lines: v m = \ ps" 1 and T p h =20 K. Dotted lines: u m =0.5 ps -1 and 
T p h=4.2K. 



case is transformed into the system 

% = -V m Xe+Ve+hu> -e B ,eXe+hu> -e B i < £ < £ B 

In the case B, when hu>o — e B < £ B < Tiuo see Fig. lc, Eq. (17) is transformed into the system 



3§ = -v m X s , e B < e < hLu Q - e B (A2) 



Joc(/|e) + Mfe+hu -e B ~ fe) + ^ e+ fc Wo _ es , e /e+ftu, -e B ~ ^ E - huo+EBt J e = 0, TlMo - £ B < £ < S B (A3) 

Joc(/|e) + vt(h-e B - 2/ e ) - v\_ huJa+eB J £ = 0, e B < e < ?iwo 



while Eq. (22) takes the form 

% = -JWe + ? E +too-M,eX E +to - eji ) < £ < TlUJo - £ B 

-§ = -(ym + vl-hu> +e B ,e)Xe + ^e+hu> -e B ,eX e+ hu l0 -e B , hu} -£ B <£ < £ B (A4) 

% = -(% + i' E -fi U|) + £ B,e)x £l e B < e < Tiuj 
Similarly for high-biased SL (case C, when hcoo < e B see Fig. Id) one obtains the system 

Jac{f\e) + ul +2nuJo _ eBi J e+ 2n Uo -e B - vl-Hu, +e B ,efe = °) < S < S B -huj 

Jac{f\e) + V t e+h ^ eB Je+hu J0 -e B - ^ e -ftu, 0+efl , e /e = °> E B - < £ < ftw 

and x E i s determined by the system 

% = -{v m + Ve-hv +eB,e)Xe) < £ < £ B - fkJ 



(A5) 



(A6) 
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